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^\ ' Abstract. We analyze the spectral lags of a sample of bright gamma-ray burst pulses observed by CGRO BATSE 

and compare these with the results of high-resolution spectroscopical investigations. We find that pulses with hard 
spectra have the largest lags, and that there is a similar, but weaker correlation between hardness-intensity corre- 
^ ' lation index, t], and lag. We also find that the lags differ considerably between pulses within a burst. Furthermore, 

the peak energy mainly decreases with increasing lag. Assuming a lag- luminosity relation as suggested by Norris 
et al., there will thus be a positive luminosity-peak-energy correlation. We also find that the hardness ratio, of 
the total flux in two channels, only weakly correlates with the spectral evolution parameters. These results are 
consistent with those found in the analytical and numerical analysis in Paper I. Finally, we find that for these 
' bursts, dominated by a single pulse, there is a correlation between the observed energy-fiux, F, and the inverse of 

, the lag. At: F oc At~^ . We interpret this fiux-lag relation found as a consequence of the lag-luminosity relation 

and that these bursts have to be relatively narrowly distributed in z. However, they still have to, mainly, lie beyond 
z ~ 0.01, since they do not coincide with the local super-cluster of galaxies. We discuss the observed correlations 
within the coUapsar model, in which the coUimation of the outfiow varies. Both the thermal photospheric emission 
Q ' as well as non-thermal, optically-thin synchrotron emission should be important. 
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1. Introduction I. We found that the spectral lag correlates strongest with 
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parameter $0 describes the decay rate or the peak energy 

analytical and numeric al simulation analysis from Paper spectrum (see eq. 2, in Paper I). For different pulses 
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the same, and thus the lag is not expected to be the same 
during the whole bursts. 

In this paper we will analyze GRB pulses, observed 
by the Burst and Transient Source Experiment (BATSE) 
on the Compton Gamma-Ray Observatory (CGRO), and 
compare the results with the trends and correlations sum- 
marized above. This will be done in Sect.|21 In particular, 
we will study the spectral lags and hardness ratios. In 
Sect. 01 we will also discuss how our results relate to the 
physical models that have been proposed to explain the 
observed lag correlations. We conclude our discussion in 
Sect.El 

2. CGRO BATSE Observations 

We defined a sample of strong and we ll-separated pulses 
by using the lBorgonovo fc R vdel ( 20(1 sample which con- 
sists of 47 pulses that all show a clear spectral evolution 
defined by the hardness-intensity correlation (HIC) dur- 
ing the pulse decay. From these we chose the cases which 
have a 'clean' rise phase that is not contaminated greatly 
by previous pulses. We also limited our sample to cases 
for which the rise time, tr, is much shorter than the de- 
cay time, td'- tr/td < 0.2. This gave us the sample of 17 
pulses in 15 bursts, which are presented in Tabled We an- 
alyzed the four-channel discriminator rates of BATSE and 
as customary, used channel 1 ('^ 25 — 50 keV) and chan- 
nel 3 (~ 100 — 300 keV). These data are constructed from 
three data types from the Lar ge Area Detecto r (LAD): 
DISCLA, PREB, and DISCSC fFishmanl. ll989l) . We also 
analyzed the 128-channel data (high energy-resolution 
data, HERB), to be able to compare the results. We decon- 
volve the observed count data into photon fluxes and per- 
formed the analysis with RMFI T 1.0b6, provided by the 
BATSE team at MSEC (see e.g. iRvde fc SvenssonI l)2000|) 
for details). 



2.1. Analysis of Pulses 

In Paper I, we derived an analytical function for a 
GRB pulse, with four parameters (see also lKocevski et all 



Table 1. Sample of 17 pulses in 15 GRBs defined by their 
peak times 



m = 



-t/tr 



(1) 



{l + {t + t,e-t/'')/TY 

Here r is the decay time scale and is connected to the 
rise phase. The decay index d is defined by requiring that 
F{t) — *■ t~'^ as t — > oo and Aq is an analytical function of 
the other parameters. We will also use the peak energy, Eq, 
at the beginning of the pulse, which is given when t = m 
Eq. Q and defines the onset of pulse emission according 
to the analytical model. The channel light-curves in our 
sample were fitted with this function. Figure ^ depicts a 
few examples of the studied bursts and their channels 1 
and 3 light curves and the fits are shown by solid lines. 
Even though the function was deduced for the bolometric 
light curve, it is flexible enough to fit the channel light- 
curves, which is reflected by the fact that the fits mainly 
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GRB911016 
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1.6 
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2.8/24 
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2 
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GRB920525 
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GRB920830 
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GRB921207 
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8.6 


GRB930201 


2156 
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14 


GRB941026 


3257 





3.0 


GRB950624 


3648 


3 


23/41 


GRB950818 


3765 


1 


66 


GRB951102 


3891 


2 


33 


GRB960530 


5478 


2 


2.0 


GRB960804 


5563 


4 


1.4 


GRB980125 


6581 





48 


GRB990102 


7293 


6 


3.2 



have good values of the . In column 2 of Tab. |2] the 

peak times are given for the light curves over the whole 
spectral range of BATSE (- 25 - 2000 keV). The peak 
time is given in seconds since the BATSE trigger, and 
is thus different from t in Eq. J^l by a constant offset. 
The main purpose of measuring this time is to compare it 
with the same time measured on the photon light-curves 
instead (see Sect. I2.3|l 

2.1.1. Analytic Lags, lagl'i 

For every pulse, the analytic spectral lag, = laglS, was 
found by measuring the time difference between the light 
curve peaks in channels 1 and 3. These peaks were iden- 
tified as the maxima of the function in Eq. (Q. These 
results are given in the third column in Tab. [5] We note 
that there is a large spread in values of the lag. The largest 
value is for the first pulse in GRB960524 (BATSE trigger 
# 3648) with 2.65 s. Furthermore, GRB 980125 (# 6581) 
and the second pulse in GRB911031 (# 973) are insignifi- 
cantly different from zero, which means that the pulses in 
the four spectral ch annels peak a t the same time. As noted 
bv Band 11997) and lNorris et~.. t2000.). GRB lags are, in 
general, concentrated towards such short time scales (< 
100 ms)). The case GRB 951102 (# 3891) is an example 
of a pulse with a negative lag, that is, the hardest radiation 
lags behind the softer radiation, even though the signif- 
icance of it being different from zero is low. Finally, the 
lags for the two cases in the sample which have two pulses 
GRB911031 (# 973) and GRB960524 (# 3648) show that 
the sizes of the lags are significantly different between the 
pulses. This is an important observation, since most often 
lags are measured over the whole burst and it is assumed 
that the lag is constant. Such a lag measurement there- 
fore more refiects the averaged lag behavior. Furthermore, 



Ryde et al.: GRB Spectroscopy 3 




Time(sec) 

Fig. 1. Light-curves of six of the studied pulses; (a) GRB911031 (# 973a); (b) GRB911031 (# 973b); (c) GRB 941026 
(# 3257); (d) GRB960524 (# 3648); (e) GRB 951102 (# 3891); (f) GRB 980125 (# 6581). The black curves are for 
BATSE channel 3 and the grey curves for channel 1. The smooth line is the best fit model (Eq. ^). 



iBorgonovo fc R vd3 1 2n0lh noted that the HIC index, r/, in 
multi-pulse bursts are often constant from pulse to pulse, 
while Liang & Kargatis (1996) noted that the HFC index, 
$0, is similarly, practically constant. Both these observa- 
tions are in contrast to that of the lags, which do differ 
significantly. F. Ryde & D. Kocevski (2004, in prep.) have 
made a more detailed analysis of the spectral lag in multi- 
pulsed bursts and showed that this is a general feature. 

2.1.2. Fluxes and Hardness Ratios 

We also calculated the photon flux [s~^ cm~^] at the peak 
of the light curve, i^pk, which appears in column 6 in Tab. 
O In Fig. 121 this flux is plotted versus the " analytic spec- 
tral lag" , lagl3. There is one data point that differs consid- 
erably from the power-law correlation that appears, and 
that is for the second pulse in trigger 973. Excluding this 
point, the best fit to the correlation is 



i^pk = 5.44 Ai 



-0.91±0.05 



(2) 



This relation is reminiscent of the lag-luminosity relation 
found by Norris et al. ( 2000) and the cor responding rela- 
tion for the photon flux ISalmonsonI l)2000() and will be dis- 
cussed further in Sect. 13. ll . The fact that the second pulse 
in 973 does not follow the general trend is of interest. As 
mentioned above the lags for different pulses vary within 
a burst. Therefore, the relation in Eq. Q (or correspond- 
ingly, the lag-luminosity relation) could either be valid 
for all individual pulses, each of which produces a point 
following the correlation, or be valid only for the domi- 
nant pulse in a burst. The latter alternative was promoted 



bv lHa,kki]afc GiblinI lEioi who discussed two bursts, of 
which a later emission component had a significantly dif- 
ferent lag. Also by removing the low-int ensity emission in 
calculating the CCF lags. lNorris et alJ l(2000() disregarded 
any differences in lag for the weaker pulses. The behavior 
of the second pulse in trigger 973, studied here, suggests 
such an interpretation as well, that is, it is mainly the 
dominant pulse that gives rise to the correlation. 

Finally, we calculated the hardness ratio (HR31), de- 
fined as the ratio of the total counts in BATSE channels 3 
and 1. We integrated the count light-curves in the interval 
which had a flux level above 10% of the peak value. The 
ratios are given in column 7 in Tab. |31 

All these measurement will be compared with the 
high-resolution spectroscopical analysis below fSect. lT^ . 
However, before doing this, we will investigate two issues 
that could affect the interpretation, namely the use of the 
cross correlation function as a measure of the spectral lag 
(Sect.l^l) and the use of count light curves fSect. lO|l . 



2.2. The Cross Correlation Function, CCF13 

It is often the case, for low signal-to-noise bursts, that 
a proper deconvolution cannot be made. Problems in fit- 
ting the channel light-curves with analytical functions also 
arise for weak bursts. Therefore, the cross correlation func- 
tion (CCF), measured on the count light-curves, has to be 
used to characterize the shifts in time. For BATSE ob- 
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Table 2. The peak times and lags of the pulses m the sample 



Trigger 


Counts 


Photon flux 


peak time [s] 


laglS [s] 


peak time [s] 


lagl3 [s] 


907 


1.55 ± 0.01 


1.3 ± 0.3 


2.1 ± 0.1 


3.1 ± 1.6 


973a 


2.84 ± 0.02 


0.75 ± 0.08 


3.1 ± 0.1 


0.8 ± 0.6 


973b 


24.04 ± 0.03 


.0004 ± 0.06 


24.1 ± 0.2 


0.1 ± 0.7 


999 


3.939 ± 0.003 


0.11 ± 0.01 


4.02 ± 0.02 


0.14 ± 0.11 


1625 


4.859 ± 0.005 


0.15 ± 0.02 


4.98 ± 0.03 


0.12 ± 0.16 


1883 


1.23 ± 0.02 


0.9 ± 0.1 


1.59 ± 0.09 


1.2 ± 0.4 


2083 


8.578 ± 0.004 


0.23 ± 0.03 


8.67 ± 0.03 


0.2 ± 0.2 


2156 


14.45 ± 0.01 


0.27 ± 0.05 


t A A 1 f\ rxf 

14.64 ± 0.05 


0.3 ± 0.3 


3257 


2.944 ± 0.001 


2.2 ± 0.8 


4.5 ± 0.3 


3 ± 5 


3648a 


23.5 ± 0.1 


2.7 ±0.6 


23.9 ±0.4 


1.7 ±4.3 


3648b 


40.99 ± 0.01 


0.4 ±0.1 


41.2 ±0.1 


0.7 ±0.8 


3765 


66.069 ± 0.004 


0.05 ±0.02 


66.12 ±0.02 


0.05 ±0.11 


3891 


33.198 ±0.005 


-0.02 ±0.03 


33.24 ±0.02 


-0.03 ±0.11 


5478 


2.00 ±0.05 


0.8 ±0.3 


2.3 ±0.2 


1.3 ±1.7 


5563 


1.418 ±0.004 


0.07 ±0.01 


1.47 ±0.01 


0.08 ± 0.06 


6581 


47.607 ± 0.002 


0.053 ±0.009 


47.66 ± 0.02 


0.08 ±0.04 


7293 


3.24 ±0.06 


2.0 ±0.9 


4.0 ±0.3 


3±4 



Table 3. The HRS parameters for the studied sample 



BATSE 
Trigger 


V 


<a> 


[keV] 


[cm-^l 


-Fpk 

[photon/s/cm^] 


HR 




907 


2.34 ±0.15 


0.37 ±0.08 


504 ± 75 


20 ± 1 


4.9 ±0.3 


1, 


.87 


0.40 


973a 


1.45 ±0.21 


-1.01 ±0.03 


515 ±25 


55± 11 


12.4 ±0.7 


1, 


.30 


0.35 


973b 


0.84 ±0.1 


-1.38 ±0.06 


544 ± 58 


17 ±4 


6.8 ±0.2 


0, 


.97 


1.7^ 


999 


1.9 ±0.1 


-0.5 ±0.1 


373 ± 35 


24 ± 1 


20.9 ± 1.5 


1, 


.81 


0.67 


1625 


2.37 ±0.3 


-0.74 ±0.02 


1005 ± 45 


42 ±7 


55.3 ± 1.1 


2, 


.35 


1.8 


1883^ 


2.39 ±0.29 


-0.3 - -2.0 


451 ± 21 


18 ± 1 


8.7 ±0.8 


1, 


.23 


0.45 


2083 


1.78 ± 0.04 


-0.36 ±0.02 


1010 ± 60 


65 ±2 


86.0 ± 1.5 


0, 


.55 


0.18 


2156 


1.55 ±0.15 


-0.84 ±0.02 


475 ± 65 


42 ±6 


23.1 ± 0.7 


2, 


.08 


0.41 


3257 


2.8 ±0.2 


0.07 ±0.06 


344 ± 22 


75 ±4 


3.8 ±0.4 


1, 


.65 


0.38 


3648a 


0.64 ± 0.03 


-0.61 ± 0.2 


157 ± 19 


19 ± 1 


1.5 ±0.4 


0, 


.53 


1.11^ 


3648b 


1.42 ± 0.08 


-0.36 ±0.05 


716 ± 17 


19 ±2 


8.5 ±0.4 


1, 


.28 


0.38 


3765 


2.4 ±0.2 


-0.76 ±0.02 


381 ± 17 


64 ±8 


45 ±2 


1, 


.86 


0.64 


3891 


1.8 ±0.18 


-0.75 ±0.06 


269 ± 13 


17 ± 1 


29.6 ± 1.2 


1, 


.36 


0.68 


5478 


1.87 ±0.16 


0.13 ±0.13 


813 ±32 


17 ±2 


3.8 ±0.9 


1, 


.30 


0.53 


5563 


2.35 ±0.3 


-0.87 ±0.07 


320 ± 15 


19 ±4 


39 ±2 


1, 


.35 


0.76 


6581 


1.62 ±0.2 


-0.65 ±0.03 


762 ± 40 


20 ± 1 


87.2 ± 1.5 


2, 


.04 


1.16 


7293 


2.04 ±0.12 


0.79 ±0.09 


397 ±51 


38 ±4 


3.7 ±0.5 


1, 


.97 


8.6 



Lag-estimated. 
^ Strong a-evolution. 
^ Secondary pulse; ignored. 



servations the spectral lag^ is often calculated between 
channels 1 (- 25-50 keV) and 3 (~ 100-300 keV). As the 
light curves (or time series) are sampled at a number of 
discrete data points sums are calculated instead of inte- 



^ Formally the lag is the variable of the correlation function 
and therefore a better notation of the time-shift between the 
pulses is the peak lag. 



grals (see eq. 4 in Paper I). The discrete CCF analysis 
was implemented by methods s imilar to those employed 
by Ifiand {199^ and iNorris et a l. ( 2000) and we denote 
this measured lag by CCF13. We therefore need to inves- 
tigate whether there is any difference in the measurements 
between analytic lags, that is laglS, and CCi^lS. 
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0.1 1 

Analytic Lag [s] 

Fig. 2. Peak photon flux as a function of spectral lag. 
The best fit is marked by the solid line and approximately 
corresponds to an inverse relation between flux and lag. 
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Fig. 3. Analytic spectral lag, At, versus the CCF lag 
CCF13. A good linear correspondence exists, as expected. 



We calculated the CCF(At, chl, ch3) and found its 
maximum value, to identify the spectral lag. At, between 
the light curves in channels 1 and 3. We fitted a cubic 
polynomial to the peak of the resulting discrete CCF func- 
tion. The statistical error in this measurement was esti- 
mated by using a Monte Carlo routine in which Poisson 
distributed noise was added to both channels individu- 
ally by an amount consistent with the observed signal-to- 
noise ratio. Each iteration produced an independent lag 
measurement, allowing for the accumulation of a cross- 
correlation peak-distribution (CCPD). This process was 
continued until the resulting CCPD approached a normal 
distribution, which typically took 250-500 runs. To this 
resulting distribution a Gaussian function was fit from 



Table 4. Spectral lags between channels 1 and 3 measured 
by their light curves {laglB) and with the CCF (CCF31). 



Spectral lags 



Trigger 


laglS [s] 


CCF31 [s] 


on? 


1 .0 31 U.O 


1 n? + n 1 1 




n 7^; + n ns 

\J. 1 o zn u.uo 


\J.OO m U.UU 


y t ou 


1 + n riR 

i-U 31 U.UU 


U.UU HI U.UU 


QQQ 


n 1 D + n ni 


U.UU 31 U.UU 






n 1 'H + n nfi 

u. J-O zn U.UU 




n 8Q + n 1 n 


n fi*? -1- n nfi 

u.uo 31 U.UU 




U.zo ± U.Uo 


U.zo ± U.UO 


2156 


0.27 ±0.05 


0.25 ±0.06 


3257 


2.2 ±0.8 


1.8 ±0.2 


3648a 


2.7 ±0.6 


1.0 ±0.3 


3648b 


0.4 ±0.1 


0.65 ±0.06 


3765 


0.05 ± 0.02 


0.06 ±0.06 


3891 


-0.02 ±0.03 


0.00 ±0.06 


5478 


0.8 ±0.3 


0.70 ±0.09 


5563 


0.07 ±0.01 


0.06 ±0.06 


6581 


0.053 ± 0.009 


0.06 ±0.06 


7293 


2.0 ±0.9 


1.8 ±0.2 



which the peak determines the most likely lag and the 
full-width half-max yields the 1 a confidence level. The 
results are also given in Tab.0] In Fig. O the analytic lag, 
laglS, is plotted as a function of the CCF lag, CCF13. A 
linear fit is given which shows that there is a good cor- 
respondence between the two methods. The best fit is 
/agl3 = (-0.02 ± 0.01) + (1.27 ± 0.08) x CCF13. The 
CCF lag measurements thus tend to underestimate the 
real time difference between the peaks of the light curves 
in the two energy bands. The bursts in our sample are 
simple pulse structures and therefore the correspondence 
might be be tter th an would be the case for very variable 
light curves. IWu fc Fcnimorc (2000) pointed out that it is 
not always reliable to determine lags with CCF. This is 
especially the case for multi-peaked bursts, for which both 
the HR and the CCF give average, quantitative descrip- 
tions. They also noted that the methods used to calculate 
the CCFs can affect the results considerably. For instance, 
the inclusion of time intervals when the signal is at back- 
ground will clearly affect the measurements. What we find 
here is that the CCF measurement represents well the ac- 
tual lag for the smooth pulses in our sample. Finally, the 
CCi^l3 for the first and second pulses in 973 and the sec- 
ond and third pulses in 3648 are again notably different, 
reconfirming the conclusion above. 

2.3. Detector Counts 

We now turn to the fact that the spectral lag analysis is 
made straight on the count rate data, which is the method 
commonly used in the literature. This introduces an im- 
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portant caveat since, in general, the count light-curves 
have a time-dependent relation to the photon (or energy) 
flux light-curves. The probability of detecting an incom- 
ing photon is energy dependent and it therefore becomes 
time dependent if the energy distribution of photons vary 
with time, that is, if there is a significant spectral evolu- 
tion. For instance, for the BATSE Large Area Detector 
(LAD) the probability of obtaining a count for an incom- 
ing photon is near zero at and below 20 keV and at 30 
keV it rises steeply. From 50 to 300 keV (and higher) the 
probability of a count approaches unity. In the same man- 
ner, the energy measurement attached to the count will be 
less correct, the higher energy the incoming photon has; 
the photon is less likely to leave 100 % of its energy in 
the detector. So, to be able to interpret the count data in 
a meaningful way, one therefore has to assume the spec- 
tral evolution has negligible effect. Indeed, assuming no 
spectral evolution during the pulse, there would be a con- 
stant correspondence between the count rate in BATSE 
channels 2-1-3 and photon flux 50-300 keV for each GRB. 
For a burst with a different spectrum, the constant would 
be different. However, significant spectral evolution does 
take place during pulses (see, e.g. Ryde (1999)) and there- 
fore the count-rate light-curve could result in a misleading 
interpretation. Note that all these effects are well-known 
and are included in the detector response matrices and are 
considered when the deconvolutions are made. 

We will therefore redo the analytic-lag analysis on the 
photon-flux light-curves instead, but still in the same en- 
ergy bands, to see how the peak times and the lags are 
affected. To be able to do this we need to deconvolve the 
detected count fiuxes through the detector response. The 
pulses in our sample are strong enough to allow us to 
do this without reducing the high time-resolution of the 
~ 64 ms DISCSC data. The detected spectra were decon- 
volve d by using the empirical GRB model I Band et al.L 
which consists of two power laws, smoothly joined 
together. The high-energy power-law, (3, was fixed to the 
averaged value, while the low-energy power-law, a, was 
left free to vary, whenever permitted by the quality of the 
data. The purpose of the fits is to make the deconvolution 
and not to find the exact parameter values and there- 
fore we can allow quite low signal-to- noise ratios. Finally, 
to arrive at the photon-flux light curve, we integrate the 
photon spectra over the relevant energy bands. 

The comparison is shown in Tab. |21 Columns 4 and 
5 contain the same information as the previous two, but 
now they are measured on the photon light-curves instead. 
In column 4 are given the peak times for the light-curves, 
integrated over the the four energy channels, and in col- 
umn 5 are given the analytic lags. The first conclusion 
to be drawn from the table is that the peak times are 
systematically underestimated by using the count light- 
curves compared to the deconvolved light-curves. The rel- 
ative difference, {tc — tp)/tc, is, in general, some 10 per 
cent, and in some cases much higher. Also the photon lags 
are, in general, somewhat larger and they also exhibit a 
slightly broader dispersion. However, even though there 



are noticeable differences they are not alarming, and es- 
pecially considering that larger errors are introduced by 
the deconvolution process, the conclusion is that the count 
light-curves, at least for the long pulses (> few seconds), 
do give a reasonable value for these studies. 



2.4. Comparison with High-Resolution Spectroscopy 

We performed the high-spectral resolution spectroscopy 
on the HERB d a ta fo llowing the method outlined in 
iRvde fc SvenssonI l)200Qi') . We fitted the instantaneous 
spectra with the highest time resolution permitted by the 
data. From these fits we calculated the time-averaged, 
low-energy, power-law index < a >. The instantaneous 
a parameter is approximately constant in all except four 
cases; triggers 999, 1625, 1883, and 2083. For instance, 
in trigger 1883, a is constant at around -0.3 for the first 
few seconds and thereafter decreases rapidly to approxi- 
mately -2. At the same time Ep decreases to ~ 30 keV. 
It is therefore likely that the observed variation in a is 
due to the limited energy band being used and that Ep 
moves beyond the low energy threshold. The low-energy 
slope could therefore very well actually be constant for 
this burst. From the evolution of the instantaneous fits 
we analyzed the peak-energy evolution, Ep = Ep(t), and 
the spectral correlations, the hardness- intensity correla- 
tion (HIC), F cx £'p' and the hardness-fluence correlation 
(HFC), dEp/dt = F/$o. The resulting parameters for the 
individual triggers are given in Tab. |31 and their depen- 
dence are shown in Figure^] Panel (a) suggests that there 
is a trend that harder spectra, that is larger a, give rise to 
larger lags. This is in agreement with the numerical sim- 
ulations in Paper I (their Fig. 3), in which such a trend is 
indeed found for all rj values. Panel (b) shows that there is 
a positive correlation, even though weak, between rj and 
a; the hard-spectral pulses have the steepest HIC power- 
laws. These two relations combined suggest that rj should 
be positively correlated with lag. Panel (c) shows a weak 
correlation in this direction. To be able to compare this 
result with the numerical simulation in Fig. 4 in Paper 
I, one must recognize the correlation between rj and a. 
This means that as rj increases so does a and the relevant 
panel changes from (d) down to (a). The dominant effect 
turns out to be the increase in lag due to a, which ensures 
the positive correlation between rj and lag discussed here. 
Finally, panel (d) shows a trend that, at least for the large 
lag cases, lower Eq gives the largest lags. The numerical 
simulation in Paper I (their Fig. 4d) shows that this is the 
case for very soft-spectra pulses (for all 77). Since, for the 
observed bursts above, a and rj are correlated [panel (b)], 
bursts with hard spectra (larger a) will predominantly 
have large 77- values. This means that the runs made in 
Paper I for large a (the earlier panels in Fig. 4 in Paper 
I; c, b, a in increasing order) are relevant mainly for the 
large 77- values, i.e., the upper section of those panels. This 
is indeed where the decrease in lag as a function of Eq 
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occurs. This will cause the behavior found in panel (d) of 
the figure above. 

In Paper I we also noted that the hardness ratios only 
weakly depend on the spectral evolution parameters, since 
their relations are complicated and the spectral parame- 
ters themselves have broad dispersions which will weaken 
any correlations. This should be most pronounced for low 
ii^o-values (see Fig. 5 in Paper I). In Fig. we plot the 
measured HRs and their relationship with the correspond- 
ing parameters Eq, a and 77. The correlations are indeed 
weak with only subtle trends. The panel on the left sug- 
gests a positive correlation between Eq and the HR, which 
is consistent with the numerical simulations in Paper I 
[their Figs. 5a and 7]. A higher initial peak energy in- 
creases the counts in the higher channels and thereby the 
HRs. The two following panels similarly suggest a posi- 
tive correlation between HR and the power-law slope, a 
and the HIC rj, respectively. According to the numerical 



simulations, harder spectra (steeper HICs) are expected 
to have larger HR [Figs. 5 b(a) in Paper I], even though 
the trends are, again, expected to be weak. A harder spec- 
trum decreases the relative count flux in the lower chan- 
nels, thereby increasing the HRs. The HR31-a correlation, 
combined with the a-rj correlation in Fig. is consistent 
with the HR3I-77 correlation in the panel on the right. 
Finally, in Fig. we plot the l ag dependence on 
/ Ft nt ■ The correlations discussed in lKocevski & Liimel 
1)200.11) emerge as expected from the analytical and nu- 
merical simulations in Paper I. The power-law index is 
0.92 ± 0.06, consistent with unity. 

3. Discussion 

Interestingly, equation (0) shows that the observed pho- 
ton (energy) flu x of a pulse i s inve rsely proportional to 
the spectral lag. iNorris et al.l l|200nl) similarly found that 
bursts with high peak-fluxes predominantly have small 
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Fig. 6. Dependence of spectral lag on <f>o/-Ftot- The outlier 
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lags (their Fig. 2). However, the flux lag relation we find 
is a much tight e r corr elation than compared to the plots 
in iNorris et al. ('2000V Our analyses differ in several as- 
pects. lNorris et al. (,200Q} study a general sample (includ- 
ing bursts with lower peak fluxes) and not only single 
pulses. The lag value for multi-pulsed events will be an av- 
eraged value and the HR rather catches the distribution of 
-Bg-values among the pulses in the studied burst. The gen- 
eral distribution in the flux-energy plane is furthermore 
determined by the HIC parameter rj. Using bursts with a 
known redshift. lNorris et al.l l 200(i[l translated their fluxes 
to luminosities and a stronger correlation emerged. More 
specifically, they found that the spectral lag, At, and the 
isotropically equivalent peak luminosity [or luminosity per 
steradian, that is, the power output of the burst assuming 
that the explosion occurs i sotropically ] , L/fl , [erg/s] (see 
also lSchaefer et al.l l|2n0lh : ISalmonsonl pOOdf) ) follows the 
relation 



L/n = 2.9 X lO^^erg/s ( ^ 



~1.14±0.20 



(3) 



Using the photon-numb er peak- luminosity [photons/s] in- 
stead, ISalmonsonI l|2000l) found the exponent to be —0.98, 
which all together is suggestive of inverse linear propor- 
tionality. We will in Sect. 13. f I below argue that Eq. ^ 
most probably is a consequence of this relation. 

A positive correlation between the peak energy 
of the spectrum, Ep, and the luminosity has been 
suggested empirically lEIovd-Ronning fc Ramirez-Rui j. 
I2.OO2; Am ati et all l2002|) . Combining this with the lag- 
luminosity anti-correlation (Eq. ^ an anti-correlation is 
inferred between peak-energy and lag. This is indeed 
consistent with the behavior of the pulses in our sam- 
ple, shown in Fig. especially for larger lags. The 
bursts with the longest lag have the lowest peak ener- 
gies. Furthermore, the positive correlations between a and 
7/ with lags (Figs. 0Ji and c) similarly imply that harder 
spectra (large a) as well as steeper HICs (large ry) are 
associated with bursts of lower luminosity. 

We also noted that different pulses within a burst 
have different lags. Such a conclusion was similarly drawn 
by Hakkila et al . (2004). Finally, we also note that 
iNorris et al.l l)2000|) found that bursts with large HRs pre- 
dominantly have small lags. For our pulses in Tab. El we 
similarly find that the largest HR bursts mainly have small 
lags. 

3.1. Origin of the Flux-Lag Correlation 

We find, in our sample, that the observed flux dispersion is 
~ 1.5 orders of magnitude and follows a lag- flux power-law 
relation (Eq. [2]). Comparing it w ith the similar rela tion 
between the luminosity and lag of Norris et al.' (2000") we 
can, in principle, have two different possibilities. First, the 
lag-flux relation is a consequence of the lag-luminosity cor- 
relation. In this case, the latter relation has to be valid for 
our sample as well. Second, Eq. (2) could be dominantly 
a consequence of the distribution of the 15 bursts in red- 
shift. In this case, Eq. [4] may, but need not, be correct. 
Since the dispersion in the measured fluxes is 1.5 orders of 
magnitude, it is obvious that for the former case a smaller 
dispersion in z is needed. Below, we will argue that it is 
indeed the lag-luminosity relation that is revealed. 
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Fig. 7. Sky distribution in super-galactic coordinates for 
the 15 BATSE bursts in the studied sample. These 
are characterized by having long uncontaminated FRED 
pulse(s). The dipole moment of the distribution is very 
large. 

The quantity measured by BATSE is the photon count 
rates and the derivation of the energy fluxes will by ne- 
cessity be model dependent. Therefore it is more useful to 
use the more fundamental quantities given by the photon 
rates in discussing these matters. The photon luminos- 
ity, Liso [photon/s], that corresponds to a certain photon 
flux value, F [photon/s/c m^j, observed f rom a source at 
a redshift, z, is given bv llCarroll et all Il992(l (see also 
iMeszaros fc MeszarosI l)l99flV ) 

= r , 

Ho Jo ^(1 + + ^Mz') - z'(2 + z')CIa 

where Hq is the Hubble constant, 51 7\/ is associated with 
the present day matter density, fl\ with the dark energy 
density, and ^Im + — 1, that is, we consider a flat 
universe at once. The luminosity distance is denoted by 
dh- The extra (1 -I- z) factor in the denominator in Eq.Q) 
account for the fact that we are using the photon flux 
and not the energy flux which is ass umed in the definition 
of dj (see, e.g, end of section 2 in IMeszaros fc MeszarosI 
Furthermore, in Eq. ^ we have neglected the K- 
correction, which takes into account the effect of z on the 
bandp ass. As shown by, for instance. iBloom. Frail, fc Saril 
l)200l|) . the quantitative distribution of fluxes does not 
change dramatically. 

For the flux-lag correlation to emerge, the dispersion 
of the quantity (f -|- z)/d\ should then be smaller than 
~ two orders of magnitude. If this is the case, then the 
dispersion in luminosities should dominate the dispersion 
of flux. It has actually already observationally been con- 
firmed that such dispersions exist. For the luminosities 
of bursts with known rcdshifts there is a dispersion of 
roughly t\yo-to- three orders of magnitude (e.g.. Fig. 6 of 
iMeszaro ? ("2001')). Taking typical values of U,m — 0.27 and 
r^A = 0.73 (Riess et al,.. .2004)). and if the pulses originate 
in bursts that are, say, from 0.5 < z < 1.0, then the ratio 
of minimal and maximal c^l is 5.5, and hence (1 -I- z)/d\ 
varies by a factor ~ 20. Note that for this ratio the exact 



value of -ffo is unimportant. With such a distribution in 
z, the lag- luminosity relation of Eq. [4] will be dominant, 
and the lag-flux correlation will be observed. The actual 
distribution in z is not known, but more arguments can 
be found to support such a range. 

First of all, it must be recognized that, remarkably, 
the quantity (1 -|- has the largest disper sion for the 

smallest z-values ijMeszaros fc MeszarosI Il988l) . If our an- 
alyzed bursts were at say z ~ 0.005 one would thus ex- 
pect a large dispersion diluting the lag-luminosity correla- 
tion. Furthermore, in the case of these small redshifts one 
would also expect a concentration of objects towards the 
galaxies in the local universe. These are distributed in a 
major planar structure centered around the Virg o cluster 
of galaxies, called the super-galactic p lane, SGP ('Pee blesL 
llflQ^tlLahav et alll200ntlNorriF[l2002i) . If the bursts in our 
sample were to be associated with these galaxies, their dis- 
tribution in super-galactic coordinates would have a high 
quadropole moment, and the conclusion would be that 
they lie within a redshift of z ~ 0.01. We note here that 
only one GRB has been detected at a similarly small dis- 
tance (GRB 980425/SN98bw, z ^ 0.0085). We therefore 
plot in Fig. Ul the sky distribution of our sample in super- 
galactic coordinates. We calculated the dipole, D =< 
cos{sgb) >, and quadropole, Q —< sm^{sgb) — 1/3 > 
moments of the angular distributions. For the 15 bursts, 
D = 0.362 (notably large) and Q = 0.0455. These bursts 
are avoiding the south super-galactic region, hence the 
large dipole moment. Since the belt around the super- 
galactical equator, up to sgb = ±30 deg in super-galactical 
latitudes, covers half the sky, one would expect a statis- 
tically significant over-density in this belt for objects as- 
sociated with the SGP. However, the quadropole moment 
we find is relatively low; in the sgb — ±30 deg-belt there 
are 8 GRBs, in the remaining part of sky 7 GRBs. In 
the BATSE exposure function, in the super-galactic plane, 
there will be higher, non-vanishing moments, albeit they 
will be small, and at least for 15 events the shot noise 
will cover them. We also produced 1000 Monte-Carlo cat- 
alogues with 15 events to test the robustness of this result. 
These catalogues follow the BATSE sky exposure. If we 
order the simulated distributions from 1 to 1000 in in- 
creasingly anisotropic distributions, the observed dipole 
value lies at 993, that is, the confidence limit is about 
99.3%. Our sample is thus not strongly distributed along 
the super-galactic plane and we therefore conclude that 
our sample is not completely within this distance, but be- 
yond it. In this connection it should be noted that iNorri^ 
(2002) studied a (totally different) sample of long-lagged 
burst and plotted them in super-galactic coordinates and 
did indeed find a large quadropole moment. This led him 
to suggest that his sample is associated with SNe type 
Ib/c in the local super cluster of galaxies. 

Second, very large z-values for the 15 GRBs in our 
sample can also be excluded. Our sample is limited to 
the brightest BATSE bursts (> 2 photons/s/cm^) and to 
cases with a single or a few well shaped pulses. The first 
constrain leaves bursts from the BATSE catalogue which 
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are normally totally dominated by lags which are very 
short ( ^ 250 ms; see Fig. 2 in iNorri j ()2002|) '). According 
to Eq. ||2J) these short-lag bursts correspond to lu minous 
explos ions leading to the high fluxes observed. The lNorri^ 
()2002|) sample also has the property that bursts with 
lags greater than ~3 s all have peak fluxes lower that 2 
photons/s/cm^. However, the measured lags in our sam- 
ple (Tab. 0J arc dominated by long lags and are thus 
low- luminosity bursts (accordin g to Eg. El ). Following 
the nomenclature introduced bv iNorri^ l)2002'l these thus 
constitute a fourth region in the peak-flux-lag plane, 
namely the bright and long-lag bursts. Since they are 
low-luminous bursts but still bright this means that they 
cannot be very distant. The second constraint that we 
used excludes the most variable bursts. According to the 
empiric al correlation b etween variability and luminosity 
( Rcicli art et alJ l)200l|) ') this again leaves low-luminosity 
bursts. To be able to study individual pulses we need 
bursts of low variability which are also bright. These cri- 
teria left us with this sample. 

Third, by assuming the lag-luminosity relation to be 
valid for our sample, we can calculate a luminosity for each 
pulse and, in combination with the measured flux value, 
we can arrive at an estimation of the redshift. These esti- 
mated ledshiits are given in the last column in Tab.|3 The 
redshifts are all below z ~ 1.8 except trigger 7293, which 
is derived to have z = 8.6. The peak of the distribution is 
at 2; ~ 0.7. In any case, this also supports the small dis- 
persion of (1 -f Note here the obvious discrepancy 
in redshifts for the multi-peaked events which arises due 
to the variation in lags, already pointed out in papers I 
and above. The sample studied is therefore most probably 
from a relatively limited distance range (low dispersion in 
z) at not very high redshifts, but mainly beyond z ~ 0.01. 

Another argument that supports the contention that 
the correlations in Eqs. (2) and (4) are the same comes 
from the cosmological time-dilation. Let us rewrite Eq. 
(4) in the form 

L/n = L'{t„/At,fr, 

where L' = 2.9 x 10^^ erg/s, to = 0.1 s and a ~ 1. For 
our purpose it is essential that a > 0, and the fact that 
Atrf is the rest-frame lag. Then the observed flux will be 
given by F = L' {to / At,t)" {1 + z) / dl = i'(to/AUs)"(l + 
z)^"'""/(i|, since Atobs — {1 + z)Ati-i. This means that one 
should consider the dispersion of (1 + instead of 

(1 + z)/d\. A positive a value will decrease the effect of 
the distance dispersion on the correlation. 

In summary, all this suggests that the correlation be- 
tween the observed fluxes and lags (Eq. [2]) is a residual 
of the lag- luminosity correlation (Eq. P]) that is still ap- 
parent. 

3.2. Origin of the Spectral Correlations 

Several observational facts appear to indicate that long 
GRBs originate when the iron-core of a massive star col- 
lapses to a black hole or a rapidly rotating neutron star, 



pro ducing a coUimated, relativistic jet (coUapsar model; 
see IWooslevI lll993|) '). The coUimation angle of the out- 
flow varies, which leads to different behaviors of the ob- 
served light curves. The observer is assumed to see the 
GRB approximately along the jet axis and the opening 
angle of the coUimation of the outflow, 9j, may vary and 
is responsible for the diversity in the quantity L/fl. Since 
is observed to be correlated to the varia bility of the 
7-ray light curve, V (|R.eichart et al.L I2OOII) . a. V - 9j- 
correlation should also exist. This was indeed shown to 
exist by Kobayashi, Ryde, & MacFadyen (2002), who used 
bursts for which an opening angle had been deduced. They 
also suggested an explanation of the correlations within 
the internal shock model. There are two possible smooth- 
ing effects and they both dependent on the bulk Lorentz 
factor: the location of the pair photosphere produced by 
synchrotron photons at R± ^ l/vT' and the angular time- 
scale tang ~ R/2cT^. If a collision happens below the pho- 
tosphere, the whole internal energy produced by the col- 
lision is converted to kinetic energy again via the shell 
spreading. It is only collisions that occur above the pho- 
tosphere that can be seen by an observer. For instance, a 
small F yields a large photospheric radius and the outflow 
then becomes more ordered and therefore the variability 
is lower. In the same manner, a small F yields a large iang 
and thus a less variable light curve. Therefore, the authors 
argued that the variability is determined mainly by the 
bulk F of the outflow and thus the V — 9j correlation re- 
duces to a F — 0j correl ation. Simulating the hy drodynamic 
collisions numerically, iKobavashi et al.l (^OOZ) found that 
the data are fitted the best with F cx 9^^ and thus the 
mass loading M oc i?/(c^F) c>c 9, since the total explosion 
energy has been fou nd to be approxima tely the same from 
burst to burst fe.g.. iFraU et all l|200lh ). This scenario is 
reasonable in the coUapsar model: a wide jet leads to a 
dim burst since a wider jet involves more mass at the ex- 
plosion, leading to lower speeds of the outflowing material. 
The change in L/fl therefore depends mainly on 6*]. The 
variability, V , is determined by the opening angle which 
in its turn also affects the quantity, L/fl cx Lconst/Of- In 
this model L ocT^. 

If the spectral lag is also determined by the angular 
spreading time scale, which is the key time scale in this 
model, the lag-luminosity relation can be found. In this 
model luminosity is scaled by the opening angle L oc 9~^ , 
while the angular spreading time is tang oc F"^ cx 9^. This 
leads to L cx t~^g, which is Eq. Q. According to Paper 
1, the lag is mainly determined by the pulse time scale T, 
which is defined as T = ^oEo/Fq. Since lag and luminos- 
ity are correlated, $o/flux should also be approximately 
correlated, as shown in Fig. |B1 

In the description above it is assumed that the varia- 
tions in Lorentz factors are mainly due to variations in the 
coUimation opening-angle of the outflow, and that the jets 
have a uniform profile out to the opening angle. However, 
an alternative scenario, for which the discussion above 
is also applicable, is the structured jet in which a cer- 
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tain b eam profile is assumed fe.g. lRossi. Lazzati. fc Reesi 
In such a case the most luminous part of the jet 
has a higher Lorentz factor. A large jet opening angle gives 
the possibility for large off-axis views. The probability for 
a certain viewing angle is indeed proportional to the view- 
ing angle itself, and thus large angles are favored. If there 
is a typical jet profile in all bursts, larger opening angle 
outflows are most often seen more off-axis, with a lower F 
as a cons e quence . 

iRvdfj ii2nn4l) found, by analyzing spectrally hard 
BATSE pulses, that the spectra can be modelled by 
a thermal black-body superimposed on a non-thermal, 
optically-thin synchrotron spectrum (see also Ghirlanda 
et al. 2003). The relative strengths of these two compo- 
nents vary. A large a-value is measured if the thermal 
component is dominant, while a lower value, more like a 
synchrotron value, is found if it is less prominent and the 
non-thermal component dominates. Small outflow veloci- 
ties, F, lead to larger photospheric radii, and less dissipa- 
tion of energy in shocks (for example) which enhances the 
thermal emission, thus leading to a harder a, as well as 
larger lags. So, a low F leads to low luminosity bursts (due 
to larger photospheric radii, as well as lower F) and hard 
spectra. This is indeed the trend in Fig. assuming that 
the lag can be translated to luminosity through the lag- 
luminosity relation. Here, that relation is observationally 
based on Fig. 2 and all the relations in Fig. ^ then cor- 
respond to relations with luminosity. A pure black-body 
emission should have a HIC of = 4. If the photospheric 
emission is dominant, a high j] should be correlated with 
a lower luminosity (Fig. ^). This also provides a natural 
explanation for a a — correlation in Fig. 

Further, Fig. ^ suggests that for low-luminosity 
bursts there is a positive luminosity-to-peak-encrgy cor- 
relation. For the smallest lags (highest luminosities) the 
peak-energies decrease and there is an anti-correlation. 
In the optically-thin synchrotron-shock model an anti- 
correlation arises most naturally. Ep cx TB'^^^ where B' is 
the comoving magnetic-field strength and 7e is the charac- 
teristic Lorentz factor of an emitting electron, and since, 
in the internal shock model, 7^ is mainly dependent on the 
relative Lorentz factor of the colliding shells and not on 
F, Ep is proportional to VB' . This is the Lorcntz-boosted 
magnetic field strength, of which the square should be 
proportional to the fireball energy density in the observer 
frame, U cx L/E?, where i? cx F^ is the typical collision 
radius. Here we have assumed that the magnetic energy 
density is a constant fraction of the total energy density. 
Therefore we have Ep cx L^/^i?"! cx L^/^F-^. With, for 
instance, L cx F^ from above, we have Ep cx L^^/^, that 
is, an anti-correlation. This is for synchrotron radiation. 
The peak energy of inverse-Compton scattered photons is 
at E^ = ^'^Ep^'^'^^ and th us it has the same dependence. 
IZhang &: Meszar'o^ l|2002l) also discussed the luminosity- 
peak-energy correlation for various other outflow models: 
For Poynting-flux dominated models Ep cx F, while for 
emission from the baryon photosphere, during the shell 
acceleration phase, Ep cx L^^"*. The last two outflows thus 



easily reproduce positive correlations. Furthermore. lRvd3 
(2004) motivated that the time-evolution of the temper- 
ature decay of the thermal pulses they studied can be 
explained either by a baryonic photosphere during the ac- 
celeration phase or by a Poynting-flux dominated wind. 
In both these models a positive correlation is, as men- 
tioned, expected. We ca n there f ore in terpret Fig. 0] fol- 
lowing the suggestions in iRvdd l)2004j) . The positive cor- 
relation for low luminosities is for pulses that are domi- 
nated by the thermal emission. The negative correlation 
(at small lags) arc then from pulses in which the non- 
thermal, synchrotron component dominates and thermal 
emission is less dominant. This is indeed consistent with 
Fig^ in which the small lag pulses have soft spectra, 
while the large lag cases have harder spectra. 

This simple description can thus accommodate the 
main correlations and trends presented above. A burst 
does not necessarily need to have the same lag for all 
pulses within it, since it is mainly internal properties that 
determine the characteristics. 

3.3. Alternative Models 

Throughout this paper it has been emphasized that the 
-Ep evolution has an important role in determining the 
spectral lags. The HFC, Ep — —F/^q, which governs this 
evolution, could be due to thermal processes like saturated 
Comptonization without extra heating, which is difficult 
to achieve in general however. Such a scenario has been 
discussed by L iang & Kargatis (1996) and Liang (1997). 
ISchaefed tno4) argued further that the lag-luminosity re- 
lation can be explained by this behavior. This model as- 
sumes an impulsive heating of a confined plasma which 
starts to cool. There should not be any continuous heat- 
ing, for example, dissipation of energy through shocks. 
Furthermore, there has to be a constant injection rate of 
soft photons which gives rise to the cooling of the elec- 
trons, Te. The thermal emission ensures that the emitted 
photons have approximately the same energy as the aver- 
aged electron. The change of internal energy of the cloud 
U, which gives the luminosity, is then L = U = NT — 
NEp, where N is the number of electrons. This equation 
is consistent with the HFC. Assuming, first, that dEp/dt 
is constant for a pulse ( which is not the exact case, see 
iRvde fc SvenssonI l)2002l) l second, that the peak energy 
decays within typical range, that does not differ greatly 
from burst to burst, (clearly also an approximation), and, 
third, that the peak luminosity, Lpoak is a good measure 
of the pulse luminosity, L, then the analytical finding in 
Paper I that 



dt 



At= I i^dE 
Jae dEp 

combined with the HFC, gives that the lag 



At cx -■ 



^pcak 



(6) 



(7) 
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$0 is practically constant between pulses, and in the model 
above it is a combination of N, distance, and jet open- 
ing angle. This is similar to the lag-$o correlation, dis- 
cussed bv lKocevski fc Liand |2003,). This model has sev- 
eral advantages but needs a number of approximations 
and strong assumptions to hold. 

As mentioned above the lag-luminosity relation is in- 
triguingly close to a pure proportional ity between t he lu- 
minosity and the inverse of the lag. ISalmonsonI l|200fl) 
noted that such a relation follows from simple relativis- 
tic kinematics if one assumes a common comoving time, 
which could be due to cooling or deceleration. Assuming 
first that the viewing angle (between the jet velocity and 
the line-of sight, l.o.s.) 9 is not large, the angle-dependent 
Lorentz factor can be approximated by 2r and the 
photon-number luminosity goes as = TF'^{1 + z)~^, 
where F'j^ is the comoving value. Since a typical comov- 
ing time scale, St goes as dt'{l + z)/2T, the observed lag- 
luminosity relation emerges. The dominant effect in the 
dispersion of lags is assumed to be due to the value of T 
directed towards the observer. This could be explained in 
two different ways. First, it is natural to assume that the 
Lorentz factor is different from burst to burst. Second, if 
the outflow is in the form of a jet, the flow velocity could 
very well be dependent on the angle from its axis, as men- 
tioned above. If the l.o.s. is along the jet axis a large F 
is observed, while if the l.o.s. is different from the jet axis 
the flow velocity might be significantly lower. The trans- 
formation of the time-scale and the flux will depend on 
the viewing angle, since they a re affected by the angle- 
dependent Lorentz factor, Pf^). llo'ka fc Nakamural 1 2Q01 1 
even argued that it is indeed the viewing angle that causes 
the dispersion in spectral lags and they showed that they 
could reproduce the observed relation between the peak 
lu minosity and lags. 

iDaigne fc MochkovitchI l|2003h presented a pulse de- 
scription within the internal shock model, which can re- 
produce the lag-luminosity relation. The peak energy is 
parameterized as Ep cx p^e^T, where p is the post-shock 
density and ec^ is the post-shock dissipated energy per 
unit mass. Standard synchrotron radiation, assuming clas- 
sical equipartition, has a; = 1/2 and y — 5/2. However, 
if the equipartition between the magnetic field and par- 
ticle energies vary with either p, e, or both, then other 
values of x and y can exist. The authors found that the 
observed spectral pulse behavior is best reproduced with 
X — y — 1/4. Their model also predicts that the hardness- 
ratio should decrease with spectral lag, which is similar to 
what is found bv iNorris et al.. (,2000.) . 



4. Conclusions 

We have studied the spectral lags in a sample of bright 
pulses, al l belonging to the c l ass of lon g duration bursts 
(see, e.g. iBalazs et al.] l|2003l) : iHorvathI ^9&) ). We find 
that there is a large spread in the measured lag values. In 
particular we find that the lag values differ considerably 
between pulses within a burst. This is important to note. 



since often the whole light curve is used to measure the 
spectral lag, and the lag measurement then by necessity is 
only an averaged value. Also, models which cannot accom- 
modate a variable lag over a pulse, that is, models which 
rely solely on external effects, need to be refined. 

We also find that a correlation between the observed 
flux and the spectral lag is valid for prominent pulses, very 
similar to the lag-luminosity relation that has previously 
been reported for entire bursts, that is, not only for pulses. 
We draw the conclusion that the observed relation is a 
residual of the physical lag-luminosity relation. The sam- 
ple does not show a significantly strong clustering towards 
the super-galactic plane, which indicates that the bursts 
are at a larger distance than z ^ 0.01. However, they 
cannot lie at very large distances, since they are mainly 
low-luminous and bright bursts. 

We have compared the analytic lags with the values 
found by using the CCF. The CCF underestimates the 
lags somewhat, but there is a linear relation between the 
measurements. Also we find that the results found by us- 
ing the count light-curves and the photon-flux light-curves 
do not differ significantly and thus permit their use. 

The observed correlations are similar to the general 
analytical and numerical simulation results in Paper I. 
Foremost we find that the lag correlates strongly with the 
decay time-scale of the pulse and that the relation is lin- 
ear. This means that the relationships including the lag, 
that have been identified, should be translatable to rela- 
tionships including the pulse time-scale. This time-scale 
is closely connected to the creation process and radiation 
processes in the out-flowing plasma. 

We also argue for a coUapsar model scenario, where 
the coUimation is the important factor leading to vari- 
ation in the observed parameters. The observations are 
consistent with the existence of several components in the 
spectra, thermal and non-thermal emission, varying in rel- 
ative strength. 
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